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Abstract 


We analyze the dynamics of models of warm inflation with general dissipative effects. We con¬ 
sider phenomenological terms both for the inflaton decay rate and for viscous effects within matter. 
We provide a classification of the asymptotic behavior of these models and show that the existence 
of a late-time scaling regime depends not only on an asymptotic behavior of the scalar field poten¬ 
tial, but also on an appropriate asymptotic behavior of the inflaton decay rate. There are scaling 
solutions whenever the latter evolves to become proportional to the Hubble rate of expansion re¬ 
gardless of the steepness of the scalar field exponential potential. We show from thermodynamic 
arguments that the scaling regime is associated to a power-law dependence of the matter-radiation 
temperature on the scale factor, which allows a mild variation of the temperature of the mat¬ 
ter/radiation fluid. We also show that the late time contribution of the dissipative terms alleviates 
the depletion of matter, and increases the duration of inflation. 
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I. INTRODUCTION 


There is a widespread belief that our Universe, or at least a sufficiently large part of 
it causally connected to us, experienced an early period of accelerated expansion, called 
inflation. This happened before the primordial nucleosynthesis era could take place and 
likely after the Planck period. Many inflationary scenarios have been proposed over the 


Q 


years (^see [i|j and references therein). Most of them rely on the dynamics of a self-interacting 
scalar field (the “inflaton”), whose potential overwhelms all other forms of energy during 
the relevant period. They generally share the unsatisfactory feature of driving the Universe 
to such a super-cooled state that it becomes necessary to introduce an ad hoc mechanism 
-termed “reheating”- in order to raise the temperature of the universe to levels compatible 
with primordial nucleosynthesis. Therefore this reheating phase appears as a subsequent, 
s^arate stage mainly justified by the need to recover from the extreme effects of inflation 
^ during which a rather elaborate process of multiheld parametric resonances followed by 


particle production Pl, y, takes place. 

As an alternative some authors have looked for inflationary scenarios leading the universe 
to a moderate temperature state at the end of the superluminal stage so that the reheating 
phase could be dispensed with altogether. It was advocated that this can be accomplished 
by coupling the inflaton to the matter fields in such a way that the decrease in the energy 
density of the latter during inflation is somewhat compensated by the decay of the inflaton 
into radiation and particles with mass. This would happen when the inflaton rolls down its 
potential, but keeping the combined pressure of the inflaton and radiation negative enough 
to have acceleration. This kind of scenario, known as “warm inflation” as the radiation 
neve. d.a.a«ca>.y. was dr. p.oposad by B.esa gQ. I* .ow ..s on 

solid grounds since it has been forcefully argued in a series of papers that indeed the inflaton 


can decay during the slow-roll (see, e.g. 


],y,Q 


and references therein). Besides, this 


scenario has other advantages, namely: {i) the slow-roll condition cjp' <^V{(p) can be fulhlled 
for steeper potentials, (ii) the density perturbations generated by thermal fluctuations may 


be larger than those of quantum origin HlilQ 
mechanism for baryogenesis [l^ . 


, and {iii) it may provide a very useful 


To simplify the study of the dynamics of warm inflation, previous works treated the par¬ 
ticles created in the decay of the inflaton purely as radiation, thereby ignoring the existence 
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of particles with mass in the decay fluid. Here, we will go a step beyond by taking into 
account the presence of of particles with mass as part of the decay products, and give a 
hydrodynamical description of the mixture of massless and non-massless particles by an 
overall fluid with equation of state p = (7 — l)p, where the adiabatic index 7 is bounded by 


1 < 7 < 2. 

On very general grounds, this fluid is expected to have a negative dissipative pressure, 
n, that somewhat quantihes the departure of the fluid from thermodynamical equilibrium, 
which we will consider to be small but still signihcant. This viscous pressure arises quite 
naturally via two different mechanisms, namely: (i) the inter-particle interactions jl^ . and 
(a) the decay of particles within the matter fluid |l6l |. 

A well known example of mechanism (i) of prime cosmological interest is the radiative 
fluid, a mixture of massless and non-massless particles, as it plays an essential role in the 
description of the matter-radiation decoupling in the standard cosmological model 0H 


Q. 


A sizeable viscous pressure also arises spontaneously in mixtures of different particles 
species, or of the same species but with different energies -a typical instance in laboratory 
physics is the Maxwell-Boltzmann gas [^. One may think of H as the internal “friction” 
that sets in as a consequence of the diverse cooling rates in the expanding mixture, something 
to be expected in the matter fluid originated by the decay of the inflaton. 

As for mechanism (ii), it is well known that the decay of particles within a fluid can be 
formally described by a bulk dissipative pressure H. This is only natural because the decay 
is an entropy-producing scalar phenomenon associated with the spontaneous enlargement 
of the phase space (we use the word “scalar” in the sense of irreversible thermodynamics) 
and the bulk viscous pressure is also a scalar entropy-producing agent. There is an ample 
body of literature on the cosmological applications of this analogy -see e.g. 0, 0. In 
the case of warm inflation, it is natural to expect that, at least, some species of particles 
directly produced by the decay of the inflaton will, in turn, decay into other, lighter species. 
In this connection, it has been proposed that the inflaton may first decay into a heavy boson 


X which subsequently decays in two light fermions 'i/jd 


22j |. This is an obvious source of 


entropy, and therefore it can be modelled by a dissipative bulk pressure H. 

Our purpose in this paper is to generalize the usual warm inflationary scenario by in¬ 
troducing the novel elements mentioned above, namely the decay of the scalar held into a 
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fluid of adiabatic index 7 rather than just radiation, and especially the dissipative pressure 
of this fluid, irrespective of the underlying mechanism. We will not dwell on the difficult 
question of the quantum, non-equilibrium thermodynamical problem underlying warm in- 

27, Qv-^ut 


flation ja |23|, IM |2 


261, 


that considered in several works 


tion H). 


29 


rather take a phenomenological approach similar to 


30 


31 


32 


(which can be traced back to the early 


studies of inflation yai). Instead of adopting a model building viewpoint and looking for 
the implications of specihc assumptions, we aim at identifying typical features of models 
that yield interesting asymptotic behavior. We resort to a qualitative analysis of the corre¬ 
sponding autonomous system of differential equations using the approach developed in Q 
that allows the consideration of arbitrary scalar held potentials. We will characterize the 
implications of allowing for various forms of the rate of decay of the scalar held, as well 
as for various forms for the dissipative pressure. We consider, for instance, models with 
scalar held potentials displaying an asymptotic exponential behavior. These arise naturally 
in generalized theories of gravity emerging in the low-energy limits of unihcation proposals 


such as super-gravity theories or string theories 


On the one hand, after the dimen¬ 


sional reduction to an ehective 4-dimensional space-time and the subsequent representation 


of the t 
nential 


leories in the so-called Einstein frame typical polynomial potentials become expo- 


37 


38|. On the other hand, the theories are then characterised by the existence of a 


scalar held that couples to all non-radiation helds, with the coupling depending, in general, 
on the scalar held. The simplest example of these features can be found in the so-called 
non-minimal coupling theories. We provide a classihcation of the relevant global dynamical 
features of the cosmological model associated with those possible choices. A limited account 
of some of the results of the present work was reported in 


One question we address is whether non-trivial scaling solutions 


33, 


35 


4G 


4lj | (hereafter 


simply termed scaling solutions) exist, i.e., solutions where the ratio of energies involving 
the matter huid and scalar held keep a constant ratio. Another class of solutions refered in 
the literature as having a scaling asymptotic behavior are those for which both the energy 
density of the scalar held and that of the matter huid decay with diherent power laws of the 


scale factor of the universe 




m 


46l | . In this latter case one of the components eventually 


dominates and thus the ratio of their energy densities becomes evanescent, in clear contrast 
to the case of the non-trivial scaling solutions. We shall term these solutions as trivial scaling 
solutions to contrast them with the previous ones sometimes dubbed tracker solutions. The 


5 

































trivial case arises in association with scalar field potentials of a power-law type and, as we 
shall see, they occur when the scalar field decays have the the same type of time-dependences 
as those required by the (non-trivial, tracking) scaling solutions. 

One of the reasons why non-trivial scaling solutions are important is that they provide an 
asymptotic stationary regime for the energy transfer between the scalar field and radiation. 
This stationary (sometimes termed “quasi-static”) regime is an assumption in the standard 


treatment of warm inflation 


ll| to evaluate the temperature of matter in the final stages. 


On the other hand, introducing this class of solutions in the kinetic analysis of interacting 
fluid s |42[ Idllj l leads to an alternative to the usual T S> 3H case, generalizing the example of 
Ref. where temperature of the matter (radiation) bath is nearly constant. 

We show that this class of scaling behavior depends not only on the asymptotic form of 


3 , 


the inflaton |35|1, but also on having an appropriate time-dependent rate for the scalar field 
decay. The additional consideration of bulk viscosity, besides being a natural ingredient in 
models with one or more matter components as well as in models with inter-particle decays, 
facilitates the Universe to have a late time de Sitter expansion. 

An outline of this work is as follows. Section II studies the model underlying the original 
idea of the warm inflation proposal, namely the model in which the inflaton field decays into 
matter during inflation thus avoiding the need for the post-inflationary reheating. This decay 
is characterized by a rate T which we shall initially assume to be a constant. Our results 
though will argue in favour of a varying T and we shall thus consider the case where T oc H. 
This yields late time scaling solutions whenever the scalar field potentials asymptotes to an 
exponential behavior. This happens regardless of the slope of the potential. Subsequently, 
section IIII, analyses more realistic models where a bulk viscous pressure term If is also 
present in the equation of state of matter. We first envisage the usual form fl = —3(H for 
that pressure and, subsequently, analyse a general model with both a varying rate of decay 
and a general form for the bulk viscosity fl = —3^ p" H^, where 2a + j3 = 2 on dimensional 
grounds. Finally, section IV provides a discussion of our results. 
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II. THE DYNAMICS OF WARM INFLATION 


Warm inflation with constant T 


We consider a spatially flat Friedmann-Robert son-Walker nniverse fllled with a self¬ 
interacting scalar field and a perfect flnid consisting of a mixtnre of matter and radiation, 
snch that the former decays into the latter at some constant rate F. For the time being 
we ignore the dissipative pressnre. We also neglect radiative corrections to the inflaton 
potential 3,13. The corresponding system of eqnations reads 


= p + ^ + R(0), 

0 = -(3i7 + F)0-W(0), 


( 1 ) 

( 2 ) 

( 3 ) 


where here and thronghont we nse nnits in which SttG = c = 1. The first two are Einstein’s 
eqnations, the third describes the decay of the inflaton. From these, it follows the energy 
balance for the matter flnid. 


p = —87 H p + . 


( 4 ) 


As usnal H = a/a denotes the Hnbble factor. 

To cast the corresponding antonomons system of fonr differential eqnations it is expedient 
to introdnce the set of normalized variables 



,2 Vi4>) 

F 


3H ’ 


( 5 ) 

( 6 ) 
( 7 ) 


along with the new time variable N = Ina. Thns we get 


x' = X {Q- 3(1 + r)) - hF(0) 2 /^ 
y' = {Q + W{cl>)x)y, 
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( 8 ) 

( 9 ) 





( 10 ) 

( 11 ) 


r' = rQ , 
cj)' = y/6x , 


where a prime means derivative with respect to N, and the dehnitions 



( 12 ) 


and 


<5 = ^ + 7 (1 - - y^) 



as well as p/(3if^) = 1 — were used. Equation (fTT|l was hrst considered in 


and is crucial for the consideration of general potentials V (0) besides the particular case of 
the exponential potential. The function Q dehned by Eq. (USD is related to the deceleration 


parameter q = —da/a? by Q = 1 + g. 

The special case where r = 0, naturally, corresponds to the absence of interaction between 
the scalar held and the perfect huid, and it is an invariant manifold of the dynamical 
system dHD-dllD. It is appropriate to refer here its major features in order to better appreciate 


the implications of the decay of the scalar held (see Table QD- 


We distinguish the hxed points of the system into those occurring for hnite values of 0 
and those associated with the asymptotic limit, (j) —^ oo. In the former case, i.e., for hnite 
(p = (f)^, the hxed points always require the vanishing of the kinetic energy of the scalar held 
(x = 0). They are located at the origin {x = y = 0), and at (x = 0,?/ = 1), on the frontier 
of the phase space domain x^ + = 1, which is an invariant manifold. For x = 0, p = 0, 

the potential must have a vanishing critical point at 0*, a case that cannot be dealt with 
the variables in use, but it is well-know that if 0* is a minimum at the origin, then it is 
a stable point and the scale factor evolves as a{t) cx t^d37) 'ppg hxed points on 

x'^ + y‘^ = 1 are given by x = 0, g = 1 and require that W = 0. This means that they can 
only occur in association with extrema of the potential. Their stability is dehned by the 
sign of iy'(0*), where 0* is the value of 0 where W(0) (and hence W) vanishes. When V 
has a non-vanishing minimum and, hence W > 0, the critical point is a stable node. When 
V has a maximum and, hence W < 0, we have a saddle point (an unstable hxed point). 
These hxed points correspond to the de Sitter exponential behavior and are accompanied 
by the depletion of the matter component (p = 0). 
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To study the critical points that occur at 0 —>• cx) (which we shall label 0oo), we carry 
out the regularization produced by the change of variable 0 = 1/0. Then Eq. flll|l becomes 

0' = —v^a;0^ , (14) 


and the critical points correspond to either to a: = 0, as previously seen, or 0 = 0. The 

|. If E(0) exhibits some 


0OO critical points depend on the asymptotic behavior of V (0) 
non-vanishing asymptotic value we have again x = 0, y = 1 corresponding to a cosmological 
constant and, hence, to a de Sitter late-time behavior. If 14(0) asymptotes towards the 
exponential potential, say V oc with A constant, there are several possible hxed values 

dependent on the ratio between and 7 (see, for instance, for details). There are 

unstable hxed points on the invariant manifolds bounding the phase-space domain for all 
possible choices of both W = —^3/2 A and 7 , namely: (0 a matter dominated solution at 
X = 0 and y = 0, which is a saddle and corresponds to a{t) oc (ii) two solutions 

dominated by the scalar held kinetic energy at x = ±1 and y = 0 which are either unstable 
nodes or saddles, and correspond to the stih behavior a{t) oc 0 oo (0 ~ where 

Kq is an arbitrary constant dehning the scalar held initial velocity. There is another hxed 
point on the -t- = 1 boundary representing a scalar held dominated solution, when 
< 9. This hxed point is stable when < 97 / 2 , and unstable otherwise (saddle). Thus 
for > 97/2 (i.e., A^ > 87 ), there is a stable hxed point in the interior of the phase space 
domain. This latter point corresponds to scaling behavior between the matter and scalar 
held energy-densities 1 ^. Q Q Id^, Q, Q- This attractor solution is characterized by 
a(t) oc and 0 — 0 o = 

There are also trivial scaling solutions for which oc a~^ and p oc where n > m 
are positive constants, when 


^1 


1 /( 0 ) = 


n y 
mj 


^ /6 — n 


2n 


A 


(15) 


where 


2n 

zu = -. 

n — m 


(16) 


Coming back to the model that includes the interaction and thus letting r be non¬ 
vanishing, we immediately see from Eq. m that, along the r-direction, all the points are 
singular points if and only if Q = 0. For hnite values of 0, as x = 0 at the hxed points, this 
requires once more = 1 so that the singular points are associated with p = 0, i.e., with the 
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depletion of the matter component. Moreover, as in the r = 0 case, these singular points are 
extrema of the potential They correspond to a de Sitter behavior (a(t) oc 

0 = 00 constant) and are either stable or unstable depending on the extremum being a 
minimum {W' > 0) or a maximum {W' < 0). In fact, at the singular points corresponding 
to extrema of the potential V{(j)), the eigenvalues found in the linear stability analysis are 


h?; = -37, 


3(1+ r) 

1 4- 

4v/6iy'(0o) 

2 

\ 

9(1 + r)2 


(17) 

(18) 


On the other hand, we no longer have hxed points a.t x = 0, y = 0 (unless 7 = 0 which 
corresponds to the perfect fluid being a cosmological constant). This happens because the 
system then evolves along the r-axis towards r —>■ oo, a behavior that can only be prevented 
by the existence of a positive minimum of the potential V{(j)). At 0 —>■ cx) the system does 
not exhibit scaling solutions anymore. The only hxed points allowed in this asymptotic 
limit are those associated with a non-vanishing, asymptotically hat potential, which thus 
corresponds to the de Sitter exponential behavior. 

Accelerated expansion corresponds to the region of the phase space where Q < 1, so that 


371/^ — 3(2 — 'y)x‘^ > 3'j — 2 . 


(19) 


This condition does not carry any dependence either on r or 0. Thus we may restrict our 
discussion to a {x, y) projection of the phase space. The condition (fTT?|l dehnes for 1 < 7 < 2 
the region between the upper branch of the hyperbolae 37 ?/^ — 3(2 — 7 ) 0 ;^ = (37 — 2) and 
the boundary -t-= 1 of the phase space domain (see Figure Q). The asymptotes of the 
hyperbolae are y = (2 — 7 )/ 7 , and we see that, as 7 increases, the inflationary region 

becomes progressively smaller. In fact the region shrinks vertically towards the x = 0, ?/ = 1 
point and it reduces to it in the limit case of 7 = 2 . 

In Ref. ^1 the end of inflation is given by the condition ~ and this event is 
associated with the beginning of the matter (radiation) domination. As it becomes apparent 
from the above discussion, the condition for the end of inflation, Q = 1 , is more general and 
does not strictly require matter domination. Taylor and Berera’s condition 11| corresponds 
to the end of slow-roll inflation (i.e., 0 ^ ~ 0 ~ x) and is extended, in the present study, to 
general 7 -fluids as 


Pm ^ 


37-2 


Ptf) 


( 20 ) 
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The independence on r of the size of inflationary region should not though be understood 
as the interaction having no effect on inflation. From Eqs (ED and m we see that the 
eigenvalues of the linearized system at the fixed points carry a dependence on r which is 
such that it renders the minima of the potential more stable and the maxima less unstable 
(as if the potential became shallower). Thus the transfer of energy from the scalar field to the 
perfect fluid favors inflation in that the system spends a longer time in the neighborhood 
of the extrema of the potential. This is exactly what is meant to happen in the warm 
inflation scenario where it is assumed that slow roll holds and argued that r allows for 
steeper potentials than those required in its absence. As discussed in [ll|, it is a simple 
matter to see that the slow-roll condition on 0 

V' v 

~3H{l + r) ~ 

is easier to satisfy if the scalar field decays, that is, if r > 1 and much easier if r S> 1. 

The fact that r increases indefinitely in the present model is a consequence of its definition, 
and merely translates the fact that, unless the system is trapped at a non-vanishing minimum 
of ld(0), H decreases towards zero. Since this is a direct result of assuming a constant T, we 
consider in the next section a more appropriate model where T decreases as the Universe’s 
expansion proceeds. 

B. Warm inflation with T (x H 


We assume that T,^ = 'ST^H where T* is a dimensionless, positive constant. As H is 
expected to be a non-increasing function of time in an expanding universe, this is a simple 
choice for the time dependence of T^ such that the decays have a stage of maximum intensity 
(when inflation occurs) followed by a progressive attenuation until it vanishes altogether. 

Now r = T* is a constant parameter and the dynamical system reduces to the three 
equations 

x' = x[Q-Z{l + r)]-W{(t))y\ (22) 

y' = [Q + W{ct>)x]y, (23) 

(f)' = a/6x , (24) 

where Q is still given by Eq. ED We see that these equations are analogous to those of 
the r = 0 case of the previous section with a different coefficient on the linear term in x 
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of Eq. (j22I)- Thus the basic qualitative dynamical features remain the same as those found 
for that model (see Table HH). The decay of the scalar held though introduces two major 
consequences worthing to be emphasized. 

Besides the fact that the origin x = 0, y = 0 is again a hxed point associated with the 
vanishing of the scalar held’s energy and, hence, corresponds to the matter domination, the 
interaction given by a non-vanishing r has the relevant effect (already found in the constant 
r model) that the stability of the minima is reinforced and that the maxima become less 
unstable. Moreover, the scalar held decay prevents the existence of the hxed points at 
X = ± 1 , 1 / = 0 , that would correspond to a behavior completely dominated by the scalar 
held’s kinetic energy (and which was, therefore, associated with a stih behavior in the r = 0 
case). 

The other major ehect of the interaction arises when we look for hxed points with x'^+y'^ < 
1 at 0 —> cx). Now, we hnd that there are always attracting scaling solutions for potentials 
that have an asymptotic exponential behavior, that is, for potentials for which W const 
when 0 —> cx) Moreover, this happens independently of the steepness of the late time 

exponential behavior which is a remarkable ehect of the present model for the transfer of 
energy from the scalar held to the matter. 

Indeed the latter solutions are given by the roots of the system of equations 


(m — 1)(m 


a 


b 


— ru = 0 , 


(25) 


and 


where u 
dehned 


cos^ 9 = 




6(1 -|- r)‘ 




and 6 are polar coordinates, x = ^ cos 6 and y 


(26) 


^ sin 6 *, and where we have 


7 


6 (l + r )2 ’ 


(27) 

(28) 


as well as Woo = —y3/2 X. It is a simple matter to conclude that the ehect of a non-vanishing 
r is such that equations (P3|) and P)|) always have one non-vanishing root within the range of 
allowed values for ^ and for cos 6, and hence there are scaling solutions regardless of the ratio 
between and 87 . Furthermore linear stability analysis shows that the scaling solutions 
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are stable. It is important nevertheless to remark that although scaling solutions emerge 
for any ratio of A^/y, the way that the 7 e// index associated with the effective equation of 
state inducing the power-law scaling behavior is shifted from the corresponding 7 value of 
the scaling solutions in the absence of decays depend on being larger or smaller than 7 . 

Assuming the potential to be asymptotically given by lA oc the latter solutions are 
a{t) oc 0 — 00 = Int^^/"^, where A is given in implicit form by 

Notice that we can dehne 7 e// = 2/{3A). A linear expansion about r = 0 in the neighborhood 
of the scaling solution (for A^ 7 ^ 87 ) yields 


A = 


37 


1 + 


A2 


{A - A) 
V37 J 


(30) 


when A^ 7 ^ 87 . So the decays have the effect of increasing (resp. decreasing) the scale factor 
rate of expansion with regard to the r = 0 case if A^ > 87 (resp. A^ < 87 ). In particular 
we can see that the scaling behavior can be inflationary, for cases where this would not 
happen in the absence of decays. For instance, taking 7 = 4/8 and A^ > 4, the condition 
for the scaling solution to be inflationary is 1 -|- r > A^/4 > 1. Thus, in this model, the 
solutions yield endless power-law inflation even for a modest scalar held decay, provided that 
the asymptotic behavior of the potential is steep enough, i.e., A^ > 87 (> 4 in the present 
example). 

Naturally, besides the scaling solutions, there can also be hxed points corresponding to 
de Sitter behavior x = 0, y = 1, whenever the scalar held potential exhibits an asymptotic, 
non-vanishing constant value. However, when the potential is asymptotically exponential, 
that there are no hxed points on the boundary x^-hy^ = 1 at 0 oo in contrast to what happens 
in the r = 0 case. 

From a thermodynamical viewpoint, the above scaling solutions are particularly interest¬ 


ing. In a universe with two components, it can be shown j42L l48l | that the temperature of 
each of the components satishes the equation 


Ti _ d 
Ti a 


_ ^ 


+ 


TliSi 


(31) 


3H) dpi dpi/m ’ 

where i = 1 or 2 , n* denotes the number density of particles of the i-species, F* their rate 
of decay, and Tj the temperature of this component. In the important case of particle 
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production with constant entropy per particle, s* = 0, we also have (pi +pi)ri = — (p 2 + 
P 2 )r 2 . Thus, taking the hrst component to be the matter/radiation fluid and the second to 
be the inflaton scalar held, we have 

T^ = Tm/r OC Tm/r ■ (32) 

As (p + = 7(1 — — y‘^)/2x^ is a constant in the scaling solutions, T^ = 3rif, with 

r a constant, implies Tm/r = Bahf, where a is another constant that depends both on r and 
on the location of the scaling solution. This yields a temperature of the matter/radiation 
component evolving as a power-law T oc Thus for a close to 1, the temperature 

of the matter/radiation remains quasi-static, whereas for cr > 1 (resp. a < 1), it increases 
(resp. decreases). Notice also that for cr ~ 0, we recover the temperature law for perfect 
fluids without dissipative effects. Provided we guarantee enough inflation, r need not be 
very large (contrary to what is usually assumed to facilitate slow-rolling). Indeed, the 
temperature of the radiation at the end of inflation is 

Tend = Tbeginning 6 , (33) 


where N is the number of e-foldings. Thus, a value of sigma lower but sufficiently close to 
1 has the potential to avoid a serious decrease of the temperature of the universe. As 

3a;? 


a = 


■r , 


(34) 


-xl-yl) 

at the scaling solutions, we see that r need not be very large to ensure that a ~ 1. 

Trivial scaling solutions generalizing those given Eq. (USD in the r = 0 case, arise in these 
models for scalar field potentials of the form 


1 /( 0 ) = 


^ n \ 2 /6(1 -I- r) — n' 


m/ 


2n 


A 


(36) 


where tu is still given by Eq. (USD and A = Ar=o/ Vl + r. For the potential to be positive 
one also requires 0 < m < u < 6(1 -|- r). The only difference with respect to the r = 0 
case lies in the dependence on r of the constant factor multitplying 0™, and translates the 
fact that there is now a different distribution of the scalar field energy density between its 
kinetic and potential parts. Indeed, we have that 


VW 


6(1+ r) 
n 



(36) 
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which shows that the extra damping of the kinetic energy part when r 7 ^ 0 , as expected. 
As in the r = 0 case, the possible emergence of the trivial scaling behavior in association 
with monomial potentials (that might or not be part of double wells) happens when the 
matter fluid is already dominating and is thus of a lesser importance in the context of warm 
inflation. 


III. WARM INFLATION WITH BULK VISCOSITY 

One of the main purposes of the present work is to assess the implications for warm 
inflation of the presence of a viscous pressure, If, in the matter component, so that the total 
fluid pressure is p = (7 — l)p + If. We may assume the expression If = —3(H which, albeit 
some causality caveats, is the simplest one may think of and has been widely considered in 


the literature 


51 


52,0- 


If the mixture of massive particles and radiation is taken as a 


radiative fluid, it is admissible to adopt ( oc where r denotes the relaxation time of the 
dissipative process. For the hydrodynamic approach to apply the condition tcoiH < 1 should 
be fulhlled. Since r oc tcoi (a reasonable assumption) and the most obvious time parameter 
in this description is H~^, one concludes that 11 ~ —/? p^, with 0 < /? < 1 . 

This modihes the held equations m and iD which now read 

+ 7P + n 


H = 


n 


p = -3 1^7 + -j FTp + Fc/) , 

while equations m and © remain in place and thus the dynamical system is now 


(37) 

(38) 


x' = X [Q -3(1 + r)] - hF(0)p^ 


y = 


X = 


r = 


[Q + W{<P)x]y, 

+2Q 


X 
rQ 

V6x , 


n 


(39) 

(40) 

(41) 

(42) 

(43) 


where 




(44) 
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and 


Q 


+ 7(1 - x ^ -y‘) + x 


(45) 


Apart from raising the order of the system, the main difference with regard to the previons 
cases lies in the modihcation introdnced in Q. In fact the bnlk viscosity term contribntes 
an additional term to it, and this changes the properties of some of the hxed points of the 
warm inflation model. These implications natnrally depend on the fnnctional form of fl, 
and next we consider some specihc choices. 


A. Warm inflation with bnlk viscosity If = —?>C,H 

Onr hrst choice is the “classical” assnmption already mentioned, fl = —3(H, where ( is 
a positive constant ensnring that the second law of thermodynamics holds. We also assnme 
that is constant as in Subsection III A I 

Given the dehnition of y in the dynamical system (jSSHlSI), we see that y = —(/H oc r. 
Therefore, dehning the constant f = —T/SC so that we have r = fy, not only Eq. (I41|l 
considerably simplihes, bnt also we do not need the r eqnation The resnlting dynamical 
system is 


X [Q- 3(1 Try)] -W{(j))y^, 

(46) 

[Q + W(0)x]2/, 

(47) 

xQ, 

(48) 

\/6x . 

(49) 


We see from Eq. fl48|l that this system has hxed points with vanishing bnlk viscosity, 
X = 0 , which send ns back to the cases already stndied in section [TTl The novel sitnations, 
however, arise when the hxed points occnr with y 7 ^ 0, which reqnires Q = 0 (meaning that 
H = is constant). 

For hnite valnes of 0, say at 0*, the hxed points are dehned by x* = 0, W(0*) yl = 0. So 
we have a hxed point at x* = 0 , ?/* = 0 , y = y*, 0 = 0 *, where y* is given by 

W = -7 ■ (50) 

Being associated with the vanishing of both 0 and V (0), this is, remarkably, a matter domi¬ 
nated de Sitter solntion. It is stable if the potential has a vanishing minimnm. Alternatively, 
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we have a line of fixed points given by x* = 0 and by 7(1 — yl) = —x*, in accordance to the 
Q* = 0 condition. This solution corresponds again to a de Sitter exponential behavior and 
has the remarkable feature that the energy densities of the scalar held and of the matter 
remain in a hxed proportion. It arises in association with an extremum of the potential V (0) 
and its stability depends on the extremum being a maximum or a minimum, the maximum 
being unstable and the minimum stable. The presence of r = fy will again contribute to 
render the minima more stable and the maxima less unstable, as already found in the study 
of the cases devoid of viscous pressure. 

It is appropriate to emphasize that, for these classes of hxed points, the matter energy 
density is not depleted by the inhationary behavior. This is due, of course, to the well- 
known fact that the bulk viscosity contributes a negative pressure and induces inhationary 
behavior. 

Now, regarding the hxed points arising at 0 —> 00 , the situation is similar to that consid¬ 
ered in Subsection msi In fact the functional forms adopted both by T and If prevent the 
existence of scaling solutions. The only hxed points associated with the asymptotic behavior 
of V (0) are analogous to those at hnite <p. If V (0) has a vanishing asymptotic value, we 
have again both the x = 0 and y = 0, x = de Sitter solutions dominated by matter, 
and when K( 0 ) exhibits a non-vanishing asymptotic value at inhnity, we have the x = 0 , 
X = — 7(1 — 2 /^) de Sitter solutions characterized by a constant ratio between the energies of 
matter and of the scalar held. Incidentally, one may remark that now the parameter r that 
represents the decay of the scalar held takes a hxed value determined by the r = fy relation 
at the hxed points. 

We conclude this Section commenting that if we were to assume the type of decay consid¬ 
ered in Subsection IIIBl i.e., a constant r, the features of the dynamical system would be the 
same as in the case just considered. This means that the remarkable ehect we found there 
that attractor scaling solutions would always emerge when the potential has an asymptotic 
exponential behavior is destroyed by the addition of bulk viscosity of the type 11 = —3(H. 

B. Warm inflation with general bnlk viscosity and decay terms 

We now extend our previous analysis to ascertain the implications of more general func¬ 
tional dependences of both T^ and 11 . 
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We assume, quite generally, = r(0)iF'^ and 11 = where 5 > 0, C, a 

and (3 are constants and moreover 2a + f3 — 2 = 0 on dimensional grounds. The latter 
condition on the parameters a and (3 implies that the dimensionless variable y becomes 
X = — 3“C (p/3iF^)“ and, hence, it reduces to y = —3“C (1 — — 2/^)“, since, as previously, 

we still have p/3H‘^ = l — — The 5 free parameter controls how fast T^ decreases with 

decreasing H during slow-roll inflation (see Eq. dSZD). Moreover, for S > 1, r decreases with 
decreasing H, and that for 5 < 1 it increases. 

The dynamical system is 4-dimensional and reads 


x' 

y' 

r' 

0 ' 


o: (Q-3(l + r))-W(0)|/2, 


(g + W(0)a;)2/, 

'dj:' 




X + Q{1 - 5) 


\/6 


X 


where Q is now given by 


Q 


2x^ -f 7 (1 — — 3“C (1 — 


(51) 

(52) 

(53) 

(54) 


(55) 


In what follows it seems reasonable to further assume that 0 < a < 1 so that /3 > 0 which 
amounts to having a bulk viscosity pressure whose importance diminishes with the expansion 
and with the dilution of matter. 

Inspection of Eqs. dSIHSSl) shows that it becomes possible to avoid the restrictive Q = 0 
condition previously found in the T = constant and 11 = —3^ H models that implied a de 
Sitter behavior at the hxed points. We have now a wider range of possibilities (our results 
for the scaling solutions are summarised in Table uni). 

At hnite values of 0 we hnd a line of hxed points x = y = r = 0, and another line of 
hxed points characterized by x = 0, 0 = (po, (1 — y‘^)^~°‘ = 3"(C/7 and any value of r. In the 
latter case 0o is the value of (p at an extremum of V (0), i.e., where W{(po) = 0, and in order 
to guarantee that y"^ < 1 we require that 


3“C < 7 • (56) 

(Notice that this amounts to having p -|- p -|- 11 > 0, hence ensuring p < 0, regardless of the 
ratio p/(3iF^). It is, thus, a condition akin to the usual weak energy condition). 
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The linear stability analysis shows that the singular points x = y = r = Q corresponding 
to matter domination are unstable. In fact the eigenvalues are 


— 0 

A, = I (1 - A) (7 - 3“C) 

A* = |(7-3“C)>0 

A. = 3fly^-(l + r)'|. 


(67) 

(58) 

(59) 

(60) 


and we see that \y is positive. 

Regarding the line of singular points with r 7 ^ 0, linear stability analysis reveals that, 
besides the vanishing eigenvalue associated with r (A^ = 0 ), the sability is once more deter¬ 
mined by the nature of the extremum of V{(j)). Indeed, the eigenvalue corresponding to y 

is 

\y = 2'^ yl {a — 1) < 0 , when a < 1 , (61) 

where y^ is a solution of (1 — ylY~°' = and the eigenvalues along x and 0 are given 

by 

Kt - ± 1 v^9(l + ry - i^ylW'(4,a ). (62) 

We see from the latter equation that > 0 requires that IV'(0*) < 0, that is a maximum 
at V (0*). Otherwise, in the case of a minimum of V (0), we have either a stable node (when 
0 < IT'(0o) < 9(1 -I- r)^/ (4^/6y^)) or a stable sink (when IT'(0o) > 9(1 +r)^/(4\^y^) > 0). 

As is well-known, some authors have resorted to the cooperative action of many scalar 
fields -the so-called“c inflation”-, both in cool inflation Q and in warm inflation , to get 
a sufficiently flat effective potential capable of driving power-law accelerated expansion. Here 
we note that this can be achieved with just a single inflationary field provided the dissipative 
bulk viscosity is not ignored. Indeed, this can be seen from from Eq. (P|). An increase in the 
parameter r as well as a decrease of y^ induce an effective reduction of the steepness of the 
potential at the maxima and a greater stability of the minima. Since y^ = 1 — 
a decrease in y^ translates an increase in C, within the admissible range ( 3“((’/7 < 1). The net 
effect is that the system spends a longer time in the neighborhood of a fixed point associated 
with a maximum of the potential (alternatively, the minima become more stable). This is 
helpful for setting the conditions for slow-roll inflation. In fact, this alleviates the need for 
a large rate of decay of the scalar field. We just need [3(1 -|- r)/y^^ > 1. 
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At 0 —>• cx), labelled 0oo, Eq- (El shows that we may have the usual hxed points corre¬ 
sponding to a non-vanishing, flat asymptotic asymptotic behavior of the potential (late-time 
approach to a cosmological constant) if Q = 0, F'/F = 0 and Fl'/Fl = 0 simultaneously. How¬ 
ever, we also And asymptotic scaling behavior in the case W (0) approaches an exponential 
behavior (hF(0oo) = — y^ 3/2 A, with constant A > 0), provided 





(63) 


which amounts to having, at 0 —> cx), F(0) x (E(0))~ and F must be asymptotically 
exponential. Notice that for 5 = 1 we recover the F x iF rate of decay considered in 
section dmi). The scaling solutions are then characterized in polar coordinates, x = ^ cos 6, 
y = ^ sin e, by 


cos 6** 


A 


VQ{i + ' 


(64) 


where u = is a root of the equation 


(1 — u) {a — bu) — rbu — (1 — u)" = 0 . 


(65) 


The quantities a and b were defined above. 

We see from Eqs. (IM|) and dni that there is always one (and only one) scaling solution, 
provided the condition (El holds. (Notice that this was precisely the condition that was 
required for the existence of fixed points at finite 0). Indeed, the first two terms of Eq. ()65|1 
are a second-order polynomial ^('w) with ^ 2 ( 0 ) > 0 and ^ 2 ( 1 ) < 0, so that it has one, and 
only one zero, in that interval (0,1). Thus the addition (subtraction) of the ^ (1 — u)" has 
the net effect of making the root of P 2 approach the origin u = 0, and the latter remains 
in the (0,1) interval provided (1^1 is valid. We also And that the location of the singular 
points corresponding to the scaling behavior is now closer to x‘^+y‘^ = 0, having a smaller y^o 
value than in the models without bulk viscosity. Moreover, linear stability analysis reveals 
that under the conditions (El and a < 1 the scaling solutions are stable, i.e., are attractors. 
These results mean that the late time contribution of the matter component is enhanced 
by the viscous pressure. This is a most convenient feature for the warm inflation scenario, 
since it further alleviates the depletion of matter during inflation and the subsequent need 
for reheating. 
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The power law behavior of these solutions is a oc ~ with A given in implicit 

form by 


(37^-2) A-^(l + r) =^r + 3^+“C^^“ ^-^(1 + ^) > (66) 


where r is here the asymptotic value of this parameter at the scaling solution (where, r oc 


l-i5 


). A linear expansion in both r and C, in the neighborhood of the scaling solution when 
r = 0, C = 0, and ^ 3j, yields 







(67) 


_2 _2 

37 


Naturally, these equations reduce to the Eqs. (OT?|) and (EDD when (5 = 1 and ^ = 0. As 
found in subsection ITTBl it is possible to define in the same manner a = 2/(3A). We 


see that now it might be possible to have < 37 if a is, for instance, a rational of the type 
a = m/{2n) with m(< 2n) and n integers. However for a = m/(2n + 1), A^ can be both 


larger or smaller than 37 , rspectively, yielding 7 e// smaller or larger than 7 . In other words, 
larger or smaller values of A with regard to the case without either decays or bulk viscosity. 

It is interesting to look at the modifications of the regions of the phase space that cor¬ 
respond to inflationary behavior arising from the consideration of the viscous pressure. In 
particular, it is important to assess how they depend on the choice of the parameters. From 
Figure El we see that the size of the inflationary region is larger than in the corresponding 
models without viscous pressure (models with the same 7 ). This was expected as the bulk 
viscosity term amounts to a negative pressure. We also see, in good agreement with this, 
that the size of the inflationary region decreases with decreasing as the importance of the 
viscous pressure diminishes. 

We conclude this section by briefly commenting on the trivial scaling solutions. The 
scalar held potentials that yield trivial scaling solutions when the matter huid dominates 
are still given by Fq. dSil). The requirement that r be a constant now translates into 


r((/.)oc (!/(</.)) 


S-1 

2 


( 68 ) 


and in addition we have a consistency condition 


1 / m — 7 


3 V 3C 
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(69) 













Once again, the inequalities 0<m<n<6(l + r) must be satisfied if the potential is to 
be positive. Moreover, as in the T oc H case, these scaling solutions happen when matter 
dominates so that they are not important for warm inflation. 


IV. DISCUSSION AND CONCLUSIONS 

In this work we have analyzed the dynamical implications for the warm inflation scenario 
of the existence of a viscous pressure in the matter content of the Universe. The dissipative 
pressure may arise either because the fluid in which the infiaton decays may be treated as a 
radiative fluid or because the different particles making up the fluid cool a different rates or 
because the particles in which the infiaton decays experience a subsequent decay in another 
particles species. We have adopted a phenomenological approach and have classified the 
asymptotic behavior of models associated with possible choices of the infiaton potential, as 
well as those arising from various functional dependences both of the rate of decay of the 
scalar field and of the viscous pressure on the matter/radiation component. In general terms 
we have considered T^ = r(0) H^, where 5 is a constant, and 11 = — 3 ^ p“ a < 1 

being a constant. 

Relevant asymptotic regimes arise in association with maxima and minima of the infiaton 
potential, at finite 0, and with the asymptotic exponential behavior of the potential at 0oo- 
In the latter case, we have found that the existence of scaling solutions depends on the 
form of the decay rate of the scalar field. Indeed, a necessary and sufficient condition to 
have scaling solutions is that the rate of decay T^ oc becomes proportional to H 

and, thus, r(0) is required to become asymptotically exponential, as Too oc {Voo)~ ■ In 
contrast to the scaling solutions found in the models without decays (the r = 0 models), 
here we find scaling solutions regardless of the steepness of the potential, that is for any 
combination of V/V = —A and 7. However the ratio between these values defines whether 
the effective value of the 7-index characterising the scaling behavior, and hence the behavior 
of the scalar field itself, is larger or smaller than the corresponding value of 7 for the models 
without decays. Indeed, 7e// < 7 when A^ > 87 and, conversely, 7e// > 7 when A^ < 87. 

Moreover inflation may be facilitated and is of the power-law type. On the one hand, 
inflationary behavior emerges in association with small values of the r parameter. On the 
other hand, the additional presence of bulk viscosity helps in avoiding a difficulty faced by 


22 



the warm inflation scenario that was raised by Yokoyama and Linde 


24| . Their argument 


was that if, on the one hand, to enhance slow-roll and simultaneously avoid the depletion 
of matter, one should have a sufficiently high rate of decay of the scalar field, on the other 
hand, this would make inflation stop earlier, since the transfer of energy from the scalar held 
to matter would make the conditions for the domination of the scalar held cease swiftly. 

Overall, the presence of dissipative pressure in the matter component (which arises on 
very general physical grounds) lends strength to the warm inflationary proposal. 
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FIG. 1: Inflationary region of the models without bulk viscosity. Inflation occurs in the shaded 
region between the hyperbole and the boundary of the phase space. The figure on the left depicts 
the variation of the region with 7 , whereas in the figure on the right we have taken 7 = 4/3. 


27 










y 



FIG. 2: Inflationary region of the models with bulk viscosity for 7 = 1. Inflation occurs in the 
shaded regions between the border lines and the = 1 boundary of the phase space. The 

lowest line corresponds to C = 1/3, the intermediate line to C = 1/5 and the uppest line to C = 1/8. 
We see that the size of the inflationary region decreases with Q. 





TABLE I: The properties of the asymptotic behavior of the T = 0 model. In this table, S stands 
for stable, U for unstable, MD for matter dominated, SFD for scalar field dominated, Min for 
minimum, Max for maximum, and Seal. Sol. for scaling solution. 
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TABLE II: The properties of the asymptotic behavior of the T oc H model. In this table, S stands 
for stable, U for unstable, MD for matter dominated, SFD for scalar field dominated, Min for 
minimum, Max for maximum, and Seal. Sol. for scaling solution. The 7 e// index is defined by 
7 e// = 2/{3A) where A characterizes the power-law expansion a oc 
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0 

0 

4^00 

= 0 

- 

U /MD 

0 

1 

4oo 

/,o /O 

o 

II 

S/SED 

0 < xo < 1 

0 < yo < 1 

4oc 

V ~ 

Jeff < 7 

S/Scal. Sol., A2 > 37 

0 < xo < 1 

0 < yo < 1 

4oc 

V ~ 

Jeff > J 

S/Scal. Sol., A2 < 37 
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TABLE III: The properties of the scaling, asymptotic behavior of the T oc r((/>)LL^, 11 = 
—model: In this table S stands for stable and Seal. Sol. stands for scaling solu¬ 
tion. The 7 e// index is defined by 7 e// = 2/{2>A) where A characterizes the power-law expansion 
a oc 



y* 

</> 

V{^) 

Eq. state 

Stability 

0 < xo < 1 

0 < yo < 1 

4^00 

V ~ r(</)) ~ 

7e// < 7 

S/Scal. Sol., A2 > 37 

0 < xo < 1 

0 < yo < 1 

4^00 

V ~ r((/)) ~ 

7e// > 7 

S/Scal. Sol., A2 < 37 
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